The dipole picture and the non-relativistic expansion

We study exclusive quarkonium production in the dipole picture at next-to-leading order (NLO) accuracy, using non-relativistic expansion for wavefunction. The light cone wave functions needed have typically been available only tree level, either phenomenological models or nonrelativistic limit. Here, we discuss compatibility of approach and compute NLO relativistic corrections to light-cone function gauge.